Observational information for f(T) theories and Dark Torsion 



O 

(N 

O 

Q 



O 

U 

Oh- 
i < 

o ■ 



m ■ 
> ■ 

oo : 

00 . 

rn ; 
od ' 

o 
o 



Gabriel R. Bengochea 1 'B 

1 Instituto de Astronomia y Fisica del Espacio (IAFE), 
CC 67, Sue. 28, 1428 Buenos Aires, Argentina 

In the present work we analyze and compare the information coming from different observational 
data sets in the context of a sort of f(T) theories. We perform a joint analysis with measurements 
of the most recent type la supernovae (SNe la), Baryon Acoustic Oscillation (BAO), Cosmic Mi- 
crowave Background radiation (CMB), Gamma- Ray Bursts data (GRBs) and Hubble parameter 
observations (OHD) to constraint the only new parameter these theories have. It is shown that 
when the new combined BAO/CMB parameter is used to put constraints, the result is different 
from previous works. We also show that when we include Observational Hubble Data (OHD) the 
simpler LambdaCDM model is excluded to one sigma level, leading the effective equation of state 
of these theories to be of phantom type. Also, analyzing a tension criterion for SNe la and other 
observational sets, we obtain more consistent and better suited data sets to work with these theories. 



I. INTRODUCTION 

Current cosmological observations, mainly from type 
la supernovae, show that the universe is undergoing ac- 
celerated expansion [H-Q. This accelerated expansion 
has been attributed to a dark energy component with 
negative pressure. The simplest explanation for this dark 
energy seems to be the cosmological constant. However, 
among many candidates @-[zJ> some modified gravity 
models have also been proposed based on, for example, 
f(R) theories g^H. 

Some models based on modified teleparallel gravity 
aresented as an alternative to inflationary models 
or showing a cosmological solution for the accel- 
eration of the universe by means of a sort of theories of 
modified gravity, namely / (Lt) [HJ, based on a modifica- 
tion of the Teleparallel Equivalent of General Relativity 
(TEGR) Lagrangian [23|, [24| where dark torsion is the 
responsible for the observed acceleration of the universe, 
and the field equations are always 2nd order equations. 
It was shown in [22[ that this fact makes these theories 
simpler than the dynamical equations resulting in f(R) 
theories among other advantages. Recently, in [25[ this 
sort of modified gravity theories was called f(T) theo- 
ries and some works have begun to develop in this area 

MM. 



In [36( the tension and systematics in the Gold06 SNe 
la data set have been investigated in great detail. Other 
authors, working with different SNe la sets found these 
were in tension with other SNe la sets and also with BAO 
and CMB [HHl]. In [37j], analyzing the Union data set 
0, the UnionT truncated data set was built by discard- 
ing the supernovae generating the tension by using the 
ACDM model to select the outliers. In [38|, perform- 
ing the same truncation procedure of [13] for 10 differ- 
ent models, it was suggested that the impact of different 
models would be negligible. 



In this work we present thorough observational infor- 
mation useful to work with f(T) theories by using the lat- 
est Union2 SNe la compilation released Q , the new com- 
bined parameter from Baryon Acoustic Oscillation and 
Cosmic Microwave Background radiation (BAO/CMB) 
39] (more suitable for non-standard models than the usu- 
ally used R and A parameters) , a Gamma- Ray Burst data 
set [4(| and constraints from Observational Hubble Data 
(OHD) gJQ- 

This Letter is organized as follows: in Section II we 
review the fundamental concepts about f(T) theories to, 
in Section III, analyze a criterion of tension to improve 
the study of the new data sets including BAO/CMB 
and GRBs. In Section IV we perform the truncation 
of Union2 calculating the relative deviation to the best 
fit of the f(T) prediction for each one of the 557 points 
following [23, |38[ in order to show the disappearing of ten- 
sion and establishing a new set suitable for f(T) theories. 
In Section V we add the OHD observational information 
and discuss some remarkable results and, in Section VI, 
we summarize the conclusions of this work. 



II. GENERAL CONSIDERATIONS ABOUT f(T) 
THEORIES 



Teleparallelism [23|,[24j uses as dynamical object a vier- 
bein field ei(x ti ), i — 0,1,2,3, which is an orthonormal 
basis for the tangent space at each point x M of the man- 
ifold: ej ■ ej = rjij, where rjij = diag(l,—l,—l,—l). 
Each vector can be described by its components e^, 



fx = 0,1,2,3 in a coordinate basis; i.e. 



Notice that Latin indices refer to the tangent space, 
while Greek indices label coordinates on the manifold. 
The metric tensor is obtained from the dual vierbein as 
9fiv(x) — rjij el (x) e J u (x). Differing from General Relativ- 
ity (GR), which uses the torsionless Levi-Civita connec- 
tion, Teleparallelism uses the curvatureless Weitzenbock 
connection I44|, whose non-null torsion is 
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The TEGR Lagrangian is built with the torsion ([TJ, 
and its dynamical equations for the vierbein imply the 
Einstein equations for the metric. The teleparallel La- 
grangian is m m, m , 



L T =T = S n ^ T»„ 



(2) 



where: 



,, , 5^ p T%-S;T%) (3) 



and v p is the contorsion tensor: 



p~ 2\ p~ p p ) 



(4) 



which equals the difference between Weitzenbock and 
Levi-Civita connections. 

For a flat homogeneous and isotropic Friedmann- 
Robertson- Walker universe (FRW), 



= diag(l, a(t), a(t), a(t)) 



(5) 



where a(t) is the cosmological scale factor. By replacing 
in (P), © and (0} one obtains 



-6 — = -6H Z 
a z 



(6) 



H being the Hubble parameter H — aa^ 1 . 

In these modified gravity theories, the action is built 
promoting T to a function f(T). The case f(T) = T cor- 
responds to TEGR. In an f(T) theory the spinless matter 
couples to the metric in the standard form. Therefore, 
the equations of a freely falling particle are the equa- 
tions of the geodesies. Moreover, the source in the equa- 
tions for the geometry results to be the matter energy- 
momentum tensor. In these aspects there is no difference 
with GR. If matter is distributed isotropically and ho- 
mogeneously, the metric is the FRW metric and all kine- 
matic equations (luminosity distance, angular distance, 
cosmological redshift, etc.) will be identical to the GR 
case. Any modification in the null geodesies followed 
by light rays will be exclusively in the scale factor a(t). 
Some authors have mentioned that f(T) theories are not 
invariant under local Lorentz transformations [20l |34| . 
However, if this would affect the viability of these mod- 
els is a subject which is currently being analyzed. 

The variation of the action with respect to the vierbein 
leads to the field equations, 

e-%(e S< nf(T) - e* T p uX S p ^f'(T) + 



Si ^(T)/"(T) + - < f(T) = 4 n G e/ T/ (7) 



where a prime denotes differentiation with respect to T, 
S.^ 1 " — e i p S p p ' 1 ' and T^ v is the matter energy- momentum 
tensor. 

The substitution of the vierbein ([5]) in ^} for i = = v 
yields 



Besides, the equation i = 1 = v is 

A&H 2 j" (T)H — J'{T)[12H 2 + AH] - f(T) = 16nGp (9) 

In Eqs. ([HI©, p(t) and p(t) are the total density and 
pressure respectively. 

In [22[ it was shown that when f(T) is a power law 
such as 



f(T) - T 



(-T) n 



(10) 



leads to reproduce the observed accelerated expansion of 
the universe, being a and n real constants to be deter- 
mined by observational constraints. 

From (J8)) along with (fTC R) , the modified Friedmann 
equation results to be (e.g. 



(2n + 1) a _ 8 

Qn+l H 2n ~ 3 



irGp 



(11) 



12 H z f'(T) + f(T) = 16ttG P 



(8) 



where p — Pmo (1 + zf + p ro (l + z) 4 , z is the cosmo- 
logical redshift and as it is usual, we will call fij = 
87rGpi /(3i/o) to the contributions of matter and radi- 
ation to the total energy density today. For a = the 
GR spatially flat Friedmann equation is retrieved. The 
case n = recovers the GR dynamics with cosmologi- 
cal constant. C omp ared with GR, n is the sole new free 
parameter (see [22[ for details). 

In the next sections, we will use ax 2 — XsNe + 
X 2 bao/cmb + Xgrb + Xohd statistic to find best fits 
for the free parameters Q m and n of a model given by 
(ITU1) using several data sets. The separate \ 2 °f SNe 
la, BAO/CMB, GRBs and OHD and the corresponding 
data sets used in this work are shown in Appendix B. 
In order to see whether our model is favored over the 
ACDM model, we will also use the information criterion 
known as AIC ( Akaike Information Criterion) [U |48[ . 
The AIC is defined as AIC — —2lnC max + 2k, where the 
likelihood is defined as C oc e x / 2 , the term —2 lnC max 
corresponds to the xinin an( l & ^ s ^ ne number of param- 
eters of the model. According to this criterion a model 
with the smaller AIC is considered to be the best, and 
a difference |AATC| in the range between and 2 means 
that the two models have about the same support from 
the data. For a difference between 2 and 4 this support 
is considerably less for the model with the larger AIC, 
while for a difference >10 the model with the larger AIC 
is practically irrelevant (49[. 



III. CONSTRAINING DARK TORSION WITH 
UPDATED DATA SETS 

We found interesting to analyze what would happen if 
we applied a criterion in order to study the consistency 
between data sets, a criterion more restrictive than the 
only fact that the confidence intervals overlap. To per- 
form this analysis, we adopted the criterion of consider- 
ing the existence of tension between a given data set and 
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another set constituted combining several data sets (in- 
cluding the first one) as the fact that the best fit point 
to the first data set is out of the 68.3% (ler) confidence 
level contour given by the combined data set. Similar 
criteria were adopted in their analysis by [3fr 38]. One 
could choose not to use this more restrictive criterion; 
however, we wanted to investigate its consequences of 
applying it to several data sets in the framework of f(T) 
theories. In [36j], for example, the best fits to sets and 
subsets of SNe are compared with the means of deter- 
mining if two of those are in tension or not, and how 
far from the confidence intervals lies the ACDM model. 
With our adopted criterion, we seek more physical con- 
sistency between best-fits, so the best fits do not drive to 
too different cosmological evolutions. The best fit which 
effective equation of state is of the phantom type [501 ] 
(w e ff < — 1) tells us about very different physics from 
the one that is not. Also, best fits that lie too far apart 
from ACDM model (n = or w e ff=-l) will indicate the 
need of more exotic models. 

Something important to consider is that the best fits 
to the SNe or their combination depend also on the fitter 
used to process the SNe data sets. Avoiding this type of 
tension we make sure that in most cases both best fits 
(SNe la and combined data sets) have similar results in 
the equation of state w e ft or the n parameter. In the 
wCDM (w=const) framework, for example, it has been 
shown that with the SDSSII (MLCS) data set [H HI 
both best fits suggest different cosmic evolutions while 
when relieving the tension this problem disappears |53j ] . 
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FIG. 1: Confidence intervals at 68.3% and 95.4% in the Q m — n 
plane for the UnionT observations of SNe la only (dashed 
lines) and UnionT SNe la, BAO and CMB (solid lines). We 
also show the best fit values to the observations of UnionT 
SNe la only (star), and to the combination of UnionT SNe la, 
BAO and CMB (triangle). 

Then, with this criterion, in our previous work [22[ 
existed some tension between the data from the SNe la 
Union sample and BAO [H| and CMB [Hj]. Taking 
the truncated UnionT data set from [33] we proceeded to 
evaluate the possible tension between UnionT, BAO and 
CMB. 



We found there was still tension between these data 
sets as shown in Fig. [TJ On one hand, the best fit to 
SNe la is of the phantom type and the one correspond- 
ing to the combined set is not and, on the other hand, the 
best fit to SNe la drives H in recombination to be 20% 
greater than for the best fit to SNe+BAO+CMB, and 
w e ff today to be 35% greater. Considering the men- 
tioned criterion, UnionT eliminates the tension between 
data in some models as shown in (38|, but it is not the 
case with f(T). 

Instead of continuing truncating the Union set, now 
we proceeded to use the latest SNe la data set Union2 
[3] (processed with SALT2 light-curve fitter [Hj|) and 
the combined parameter BAO/CMB. This combined 
BAO/CMB parameter implemented in [3{| is more suit- 
able to add constraints to non-standard models (see Ap- 
pendix B for details). 

The best fit to the Union2 SNe la data set only, was 
achieved with n — 0.49 and Q, m = 0.33 with the reduced 
Xmin/ V — 0-98 (or equivalently, Ayi in = -0.18 with 
respect to ACDM with Q, m = 0.27 Hf), where v is the 
number of degrees of freedom. All the results with their 
corresponding la uncertainties and the analysis from the 
AIC criterion are summarized in Table [TT1 

Working only with the new BAO / CMB parameter we 
found the value of n for the best fit is remarkably higher 
than with other data sets {n = 4.58), although more 
efficient in constraining fi m (in contrast to working with 
BAO and CMB parameters separately) having a range of 
values more consistent with other observations, such as 
weak lensing and its combination with CMB and SNe la 
(e.g. [13, HI]). In our case we found Q m = 0.28. These 
values (n = 4.58, il m = 0.28) perform a better fit than 
ACDM by a Ax^m = —1-46. Comparing these results 
with our previous work [22| , we also found that including 
this combined parameter, the value of n for the best fit 
using BAO / CMB or its combination with the rest of the 
data sets results always greater than zero. Therefore, 
the effective dark torsion is of the phantom type (Hoj . 
This result is also in opposition with recent constraints 
when BAO and CMB are used se par ately through the 
parameters A and R respectively [261 ] . Combining the 
SNe la data with BAO/CMB data we found the best fit, 
which can be seen in Table [TTJ 

We added to our analysis a data set of observations 
of Gamma- Ray Bursts (GRBs). Knowing that there are 
still debates about if these objects are standard candles 
(e.g. [1^, [6(|), we followed the policy assumed in other 
published works such as duu to see how our results are 
modified and to check the consistency of both approaches 
presented by [6l| an d libl in the framework of f(T) the- 
ories. Recently, in [621 ] it has been demonstrated that 
the data set of [i(| is consistent with Union2. One could 
do the analysis using data sets consisting of SNe la and 
GRBs separately. Otherwise, an analysis can be made 
compiling Union2 and one GRBs data sets together [6l| . 
These two separate analysis showed identical results (see 
Appendix A). 
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Here we used the approach developed in [40j firstly 
for being a set- independent data set, this means this set 
is not only applicable with Union2, but with any other 
SNe la data set. Most importantly, we used this data 
set because evaluating it separately from the SNe la is 
more helpful in our work of finding tension between SNe 
la data and other data sets. Otherwise, in case of finding 
tension we might have needed to truncate a combination 
of SNe+GRBs data set, being this last a sum of different 
data. When incorporating GRBs to the joint statistic we 
observed that the addition of the mentioned observations 
slightly reduce the size of the confidence intervals and the 
joint best fit is displayed in Table [TT1 

In Fig. [5^, we show the confidence intervals of SNe 
la with the combination of SNe+BAO/CMB+GRBs. 
There, it can be easily seen that with the adopted cri- 
terion for the tension between data sets, there is a slight 
tension between Union2 and the other data sets. 




0.26 0.28 0.30 0.32 0.34 0.36 0.38 0.40 



FIG. 2: (a) Confidence intervals at 68.3% and 95.4% in the 
fl m — n plane for the Union2 observations of SNe la only 
(dashed lines) and Union2 SNe+BAO/CMB+GRBs (solid 
lines). We also show the best fit values to the observations 
of Union2 SNe la only (star), and to the combination of 
Union2 SNe+BAO/CMB+GRBs (dot), (b) Idem (a), for the 
Union2T data set. 

In the next section, we will adopt the criterion pre- 
sented in [37] in order to perform the truncation of 
Union2 SNe la data set with the objective of dissipat- 
ing the tension between data sets when analyzing these 
f(T) type of theories. 



IV. BUILDING A NEW IMPROVED DATA SET 

We followed the simple method used in [36 l [37l| to find 
the outliers responsible for the tension. In [361 the dis- 
tance moduli of the six SNe la which are mostly responsi- 
ble for the tension in Gold06 data set differ by more than 
1.8a from the ACDM prediction. In [37| 21 SNe la were 
discarded to build the UnionT set in order to eliminate 
the tension with CMB and BAO. 



Similarly, we firstly fitted our f(T) model to the whole 
557 SNe la in the Union2 data set and found the best fit 
parameters (with the corresponding (1q = 43.15). Then, 
we calculated the relative deviation to the best fit pre- 
diction, I /i {, s — ix t h I /o'nhs, for all the 557 data points. 
We found that as in [37| the cut 1.9a solved the tension 
problem. This cut implied to take out 39 SNe la from 
Union2. With the remaining 518 SNe la we build the 
Union2T data set. The outliers are shown in Table fl] 

TABLE I: The names of the outliers from Union2 data set. 
Outliers from Union2 

1995ac, 1998dx, 1999bm, 2001v, 2002bf, 2002hd, 2002hu 
2002jy, 2003ch, 2003ic, 2006br, 2006cm, 2006cz, 2007ca 
10106, 200511, 20051p, 2005fp, 2005gs, 2005gr, 2005hv 
2005ig, 2005iu, 2005jj, 1997k, 1998ba, 03D4au, 04D3cp 
04D3oe, 03D4cx, 03Dlco, d084, el40, £308, g050 
g!20, m!38, 05Str, 2002fx 



In Fig. [2Jd, it is shown the result of using Union2T 
with the combination of Union2T+BAO /CMB+GRBs. 
The best fit to Union2T was achieved with the val- 
ues n = 0.23, O m = 0.31 and with Xmin/ 1 ' = 0-67 
(Axmin = — 0.20), whilst the best fit to the joint anal- 
ysis of Union2T+B AO / CMB+GRBs was obtained with 
n = 0.09, n m = 0.29 with Xmm /v = 0.67 (A X 2 mm = 
—0.97). From these results we can see that the Xmin as 
the Xmin/ U have been significantly improved in respect 
to the values obtained in the previous section. 

This Union2T set along with the corresponding sets 
of BAO/CMB and GRBs are more consistent between 
them, considering the adopted criterion which advan- 
tages were mentioned above. 

V. ADDING OHD DATA 

In this section we wondered about how the results of 
the previous sections were modified when a data set with 
Hubble parameter observations H(z) was added to the 
X 2 statistic. The details of the used data are displayed 
in Appendix B. 

In Fig. [3^,, the confidence intervals are shown at 68.3% 
and 95.4% for the SNe la data only and for the combina- 
tion of SNe+BAO/CMB+GRBs+OHD. Also, displaying 
the corresponding best fits. This analysis showed that the 
best fit was reached with n = 0.33 and Q m = 0.29 with a 
X 2 m inl v = 0-96 (Ax min = —2.6). It is observed that the 
confidence intervals are slightly smaller in size than when 
OHD data is not added. Surprisingly, we found that the 
adding of OHD pushes the ACDM model out of the la 
confidence level of the combined data. The inclusion of 
OHD data, then, favors an equation of state of the phan- 
tom type and makes the values that lie at 68.3% to be 
n £ [0.12,0.58], Q m £ [0.27,0.31]. A similar result was 
obtained when combining BAO/CMB, the shift param- 
eter R of CMB and supernovae la from SDSS by using 
the MLCS light-curve fitter (Fig. 1 of (H). 
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FIG. 3: (a) Confidence intervals at 68.3% and 95.4% in the 
fl m — n plane for the Union2 observations of SNe la only 
(dashed lines) and Union2 SNe+BAO /CMB+GRBs+OHD 
(solid lines). We also show the best fit values to the observa- 
tions of Union2 SNe la only (star), and to the combination of 
Union2 SNe+BAO /CMB+GRBs+OHD (dot), (b) Idem (a), 
for the Union2T observations. 



Again, with the adopted criterion about the ex- 
istence of tension between data sets we see there 
is a slight tension between Union2 and Union2 
SNe+BAO/CMB+GRBs+OHD. Performing the ana- 
logue procedure of the previous section, we found that 
a 1.9a cut was suitable to remove the tension between 
data sets. The result is displayed in Fig. [3p. 

The combination Union2T SNe la with BAO/CMB, 
GRBs and OHD allowed the ACDM model again to lie 
inside the la joint probability region. 

We summarize in Table [TT] the main results from the 
analysis performed in this work. 



TABLE II: Best fit values and la errors for each parameter 
marginalizing over the other, for the models considered in this 
work. (1) SNe, (2) BAO/CMB, (3) SNe + BAO/CMB, (4) 
SNe + BAO/CMB + GRBs, (5) SNe + BAO/CMB + GRBs 
+ OHD, (6) U2T, (7) U2T + BAO/CMB + GRBs and (8) 
U2T + BAO/CMB + GRBs + OHD. SNe stands for SNe la 
from Union2 data set, U2T stands for the truncated Union2 
data set with a 1.9cr cut and AAIC = AICf(r) — AICacdm ■ 



Data set(s) 


n 






Ay 2 ■ 


AAIC 


1 

2 
3 
4 
5 
6 
7 
8 


o.49ii:& 

4.581:/" 
0.15i°;?I 

u - iD -0.18 
qq+0.25 
U - JO -0.21 
n 23+ 108 

om±° f 8 

0.28±g : « 


n qq+0.08 
u -°°-0.19 

0.28±H 2 2 
0.29lH 2 2 

n 2q +0 02 

u - zy -0.02 

29+ ' 02 

n ,1+0.08 

u - zy -0.02 
u - zy -0.02 


0.98 

0.97 
0.96 
0.96 
0.67 
0.67 
0.67 


-0.18 
-1.46 
-0.99 
-0.93 
-2.6 
-0.20 
-0.97 
-1.83 


1.82 
0.54 
1.00 
1.07 
-0.60 
1.80 
1.03 
0.17 



VI. CONCLUSIONS 

We have updated the constraints to an f(T) = T — 
a(—T)~ n theory by using the latest type la supernovae 
data set Union2, the new combined BAO/CMB parame- 
ter, a Gamma-Ray Bursts set and Hubble Observational 
Data. 

When the new BAO/CMB parameter is used instead 
of the frequently used A and R parameters separately 
the best fit values change with respect to previous works. 
From Table im we see that all best fit values are for n > 0, 
leading the effective equation of state to be phantom like. 
Note that in all cases where SNe la data were involved, 
we used the Union2 data set which was processed with 
SALT2 fitter and this could be an additional factor in the 
obtained results as showed in [53| . Adding GRB data did 
not modify the results appreciatively and we also found 
that two approximations performed by different authors 
are consistent between them and lead to the same results. 

We found that when including information from 
OHD to put constraints, as in the BAO/CMB case, 
an equation of state of the phantom type is favored. 
Remarkably, the simpler ACDM model lies outside 
the 68.3% confidence level region of the combined 
SNe+BAO/CMB+GRBs+OHD data. The values that 
lie at 68.3% are in the ranges n € [0.12,0.58], Q, m G 
[0.27,0.31]. 

The adopted criterion of tension between data sets 
in this work and the truncation process performed to 
Union2 data set allowed us, firstly, that the physics that 
determines the cosmological evolution through the n pa- 
rameter does not differ much when only the SNe are con- 
sidered or when those are combined with other data sets, 
and secondly, that each data set is consistent amongst the 
others. In every case, eliminating tension leaded to re- 
duce at least to one half the ratio between values of w e f / 
in z = and H in recombination, obtained from the 
best fits to SNe alone or their combination with other 
data sets. Additionally, the removal of tension by this 
criterion resulted in, when combining Union2T data set 
with BAO/CMB+GRB+OHD, the ACDM model to lie 
inside the la joint probability region of the data set with 
all observations. From Table HH the values for AXmin 
always favored the f(T) models, whilst when perform- 
ing an analysis with the AIC criterion we found that 
all the evaluated cases of f(T) have the same support 
from the data with respect to the ACDM model since 
< | A ATC | < 2. 
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Appendix A: Analysis from Hymnium Gamma-Ray 
Bursts data set 



Here we show the result of adding to the performed 
analysis with the Union2 SNe la data set, the set Hym- 
nium of 59 GRBs according to (olj . 

With the objective of comparing the obtained results 
when adding GRBs to SNe la data and BAO/CMB us- 
ing the 5 values set of and when adding the data 
sets in the way performed in |6lj . we observed the dif- 
ferences in the 68.3% and 95.4% confidence intervals of 
both methods. The result using the 59 GRBs set of [6l[ 
is displayed in Fig. |4] The values in the joint best fit are 



for n = 0.15_ 



and fL 



0.29 



+0.02 
-0.02- 



As can be seen, 



the result is very similar to the one obtained in section 
III (Fig. Hi). 



parameter but it is independent of the data points. Fol- 
lowing [63], the minimization with respect to /i can be 
made by expanding the XsNe with respect to /j,q as 



where. 



^(P) 



X%Ne{v)=A-2 f i B + f jLlC 



(B3) 



JV=557 

E 

i=l 
JV=557 

E 

i=i 

N=557 



C = 



[^obs(zj) - Mt/t( z »;Mo = o>p)] 2 

lHobs(zi) - (lth(Zil = 0, p)] 



2 {Zi) 



y - J - 




Equation (|B3I) has a minimum for /io = B/C at 

B(P) 2 



XSJVe(P) = ^(P) - 



c 



(B4) 



Since XsNe,min = XsNe,min obviously, we can instead 
minimize XsNe which is independent of (j,q. 



0.22 0.24 0.26 0.28 0.30 0.32 34 0.36 0.38 



FIG. 4: Confidence intervals at 68.3% and 95.4% in the f2 m -n 
plane coming from combining SNe+GRBs+BAO / CMB data, 
where the 59 GRBs from the Hymnium data set were added 
to Union2 according to [oil ]. 



Appendix B: Cosmological constraints methods 
1. Type la Supernovae constraints 

The data points of the 557 Union2 SNe la compiled in 
Q are given in terms of the distance modulus /j, Q bs i z -i) an d 
the corresponding uncertainty for each observed value 
cr(zj). On the other hand, the theoretical distance mod- 
ulus is defined as 



fhh(zi) = 5logi D L (zi) + /i 



(Bl) 



where /io = 42.38 — blogi^h and h is the Hubble constant 
Hq in units of 100 km/s/Mpc, whereas the Hubble-free 
luminosity distance for the flat case is 



D L {z) = {\ + z) 



dz' 



E(z\p) 



(B2) 



in which E = H/Hq, and p denotes the model parame- 
ters (here, n and Cl m ). The parameter fj,Q is a nuisance 



2. Combined BAO/CMB parameter constraints 



In the f(T) theories considered here l|10p. for later 
times the term —a/(—T) n is dominant, while in early 
times when H — > oo General Relativity is recovered. 
Also, since this model presents matter domination at the 
decoupling time as the standard model, we can use the 
BAO and CMB information as showed in [22j. 

When analyzing CMB and BAO observations there 
are two parameters commonly employed, R [o4| and A 
54]. However, a more model- independent constraint can 
be achieved by multiplying the BAO measurement of 
fs(zd)/Dv{z) with the position of the first CMB power 
spectrum peak (5f| I a — ndA(z*)/r s (z*), thus cancelling 
some of the dependence on the sound horizon scale [39l ] . 
Here, ^(z*) is the comoving angular-diameter distance 
to recombination, r s is the comoving sound horizon at 
photon decoupling, Zd « 1020 is the redshift of the drag 
epoch at which the acoustic oscillations are frozen in, and 
Dy is defined as (assumed a ACDM model) [54| . 



D v (z) = 



( /„ H(z')) 



-l 1/3 



H{z) \J Q H(z') 



(B5) 



We further assume z* = 1090 from [55j (variations within 
the uncertainties about this value do not give significant 
differences in the results). 

In [65| was measured r s (zd)/ Dv{z) at two redshifts, 
z = 0.2 and z = 0.35, finding r s {z d )/D v (0.2) = 0.1905± 
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0.0061 and r s {z d )/ D v (0.35) = 0.1097 ± 0.0036. Com- 
bining this with I a gives the combined BAO/CMB con- 
straints 13911 : 



<Ia{z*) r s (z d ) 



D v (0.2) r s (z„) 
d A (z*) r s {z d ) 



£V(0.35) r.(z,) 



18.32 ±0.59 
10.55 ±0.35 



(B6) 



Before matching to cosmological models we also need 
to implement the correction for the difference between 
the sound horizon at the end of the drag epoch and the 
sound horizon at last scattering. The first is relevant for 
the BAO, the second for the CMB, and r s (z d )/r s (z^) — 
1.044 ±0.019 (using values from 55]). Inserting this into 
(|B6p and taking into account the correlation between 
these measurements using the correlation coefficient of 
0.337 calculated by (65[, gives the final constraints we 
use for the cosmology analysis [39] ]: 



A l = 
A 2 



d A {z 



D v (0.2) 

dA(z*) 



LV(0.35) 



= 17.55 ±0.65 

= 10.10 ±0.38 (B7) 



Using this BAO/CMB parameter cancels out some of 
the dependence on the sound horizon size at last scat- 
tering. This thereby removes the dependence on much 
of the complex pre-recombination physics that is needed 
to determine that horizon scale [39(. In all the cases, we 
have considered a radiation component ft r — 5xl0 -5 . 



So, for our analysis we add to the \ statistic: 



N=2 



Xbao/cmb 

(p) = E 



[A obs (zj) - Ath(zj-,p)] 

<J 2 A (Zi) 



(B8) 



where p = (n, Q m ) are the free parameters, A b s is the 
observed value (Ai and A 2 ), A t h is the predicted value by 
the model and a a is the la error of each measurement. 



3. Gamma-Ray Bursts constraints 



Following [66| and [40|, we consider the well-known 
Amati's E pi — E iso correlation [67i - [70j in GRBs, where 
E Pt i = E Pt0 i, s {\ + z) is the cosmological rest-frame spec- 
tral peak energy, and Ei SO is the isotropic energy. 

In [7l| , it was defined a set of mo del- independent dis- 
tance measurements {r P (zi)}: 



r p (zi) 



r P (z ) 



with, 



(l + z) 1 ^ H 



r{z) 



(B9) 



(B10) 



where r(z) = di,[z)/{\ + z) is the comoving distance at 
redshift z, and zq = 0.0331 is the lowest GRBs redshift. 



Following the method proposed by [7l|, in [40J were 
obtained 5 model-independent distances data points and 
their covariance matrix by using 109 GRBs via Amati's 
correlation. The resulted modcl-indcpcndcnt distances 



~data ( 



, ,, l {z.i) and their uncertainties, the correlation matrix 
and the covariance matrix are these ones from [4fj| . 

So, a given cosmological model with {p} free parame- 
ters can be constrained by GRBs via 



[Af p ( Zl )f ■ (C^ B ) • [AfpC*)] (Bll) 

(B12) 



Xgrb(p) 
Af p ( Zl ) = f d p ata {z t )-f p {z % ) 



where f p (zi) is defined by (|B9|) and Cg)j B is the inverse 
of the covariance matrix. In this way, the constraints for 
a large amount of observational GRBs data is projected 
into the relative few quantities f p ata (zi). 



4. Observational Hubble Data (OHD) constraints 

The Observational Hubble Data are based on differ- 
ential ages of the galaxies [z2]. In [73| it was obtained 
an independent estimate for the Hubble parameter using 
the method developed in (72|, and the authors used it 
to constrain the equation of state of dark energy. The 
Hubble parameter depending on the differential ages as 
a function of redshift can be written as, 



H(z) 



1 



1 



dz 
z dt 



(B13) 



So, once dz/dt is known, H(z) is obtained directly [74j . 
By using the differential ages of passively-evolving galax- 
ies, in J7J] was obtained H(z) in the range of 0.1 < z < 
1.8 and in [43| new data at 0.35 < z < 1 were studied. 
Here, Hn from [HJ and eleven observational Hubble data 
from [43| are used. 

In addition, in [4l| the authors took the BAO scale 
as a standard ruler in the radial direction, called "Peak 
Method", obtaining three more additional data (in 
km/s/Mpc): H(z = 0.24) = 79.69 ±2.32, H(z = 0.34) = 
83.8 ± 2.96 and H(z = 0.43) = 86.45 ± 3.27, which are 
model- and scale-independent. We just consider the sta- 
tistical errors. 

The best fit values of the cosmological model parame- 
ters from observational Hubble data are then determined 
by minimizing, 



iV=15 

X 2 ohd(p) = 



H bs(zi) - H th (zt;p)} 2 



(B14) 



where as before, in this work p = (n, f2 m ), H t h is the 
predicted value for the Hubble parameter, H b s is the 
observed value, cr(zi) is the standard deviation of each 
measurement, and the summation is over the full 15 val- 
ues at redshift z, mentioned above. 
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